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This paper is the first in a series of three papers that aim at understanding the scaling behaviour of 
hydrodynamic turbulence. We present in this paper a perturbative theory for the structure functions 
and the response functions of the hydrodynamic velocity field in real space and time. Starting from 
the Navier-Stokes equations (at high Reynolds number Re) we show that the standard perturbative 
expansions that suffer from infra-red divergences can be exactly resummed using the Belinicher-L’vov 
transformation. After this exact (partial) resummation it is proven that the resulting perturbation 
theory is free of divergences, both in large and in small spatial separations. The hydrodynamic 
response and the correlations have contributions that arise from mediated interactions which take 
place at some space- time coordinates. It is shown that the main contribution arises when these 
coordinates lie within a shell of a “ball of locality” that is defined and discussed. We argue that the 
real space-time formalism developed here offers a clear and intuitive understanding of every diagram 
in the theory, and of every element in the diagrams. One major consequence of this theory is that 
none of the familiar perturbative mechanisms may ruin the classical Kolmogorov (K41) scaling 
solution for the structure functions. Accordingly, corrections to the K41 solutions should be sought 
in nonperturbative effects. These effects are the subjects of papers II and III in this series, that will 
propose a mechanism for anomalous scaling in turbulence, which in particular allows multiscaling 
of the structure functions. 

PACS numbers 47.27.Gs, 47.27.Jv, 05.40.+j 


1. INTRODUCTION TO THE SERIES OF 
PAPERS 

This paper is the first in a series of three papers that 
aim for a systematic elucidation of the analytic structure 
of the the statistical theory of hydrodynamic turbulence 
at high Reynolds number Re. The natural quantities 
that appear in the statistical description of turbulence are 
averages of products of hydrodynamic fields at different 
space-time points. The aim of the analytic approach is 
to evaluate such quantities on the basis of the equations 
of fluid mechanics. 

The fundamental held in hydrodynamics is the veloc¬ 
ity held of the huid, denoted as u(r, t) where r is a point 
in d-dimensional space (usually d = 2 or 3) and t is the 
time. Statistical quantities that attracted years of exper¬ 
imental and theoretical attention p],|]|| are the structure 
functions of velocity differences, denoted as S n (R) 

S n {R) = (|u(r + R,f) - u(r,f)| n ) . (1.1) 

where (...) stands for a suitably dehned ensemble av¬ 
erage. It has been stipulated for a long time that the 
structure functions scale as a function of R when R is in 
the so called “inertial range”, i.e. 77 <C R <C L with 77 be¬ 
ing the inner viscous scale and L being the outer integral 
scale of turbulence: 


with being scaling exponents. One of the major ques¬ 
tions in fundamental turbulence research is whether these 
scaling exponents are correctly predicted by the classical 
Kolmogorov 41 theory in which ( n = n/ 3, or whether 
these exponents manifest the phenomenon of “multiscal¬ 
ing” with being a nonlinear function of n, as has been 
indicated by experiments. 

Experimental research did not confine itself to the mea¬ 
surement of the structure functions of velocity differ¬ 
ences. Gradient fields featured as well. For example the 
correlation function of the energy dissipation field has 
been studied extensively. The dissipation held e(r ,t) is 
dehned as 

e(r,f) = ^[V a U/ 3 (r,f) + Vpu a (r,t)} 2 (1.3) 

with v being the kinematic viscosity. The correlation 
function of the dissipation held K ee (R) is 

Kee(R) = (e(r + R, f)e(r, t)) , (1.4) 

where e(r, t ) = e(r, t ) — e. Here and below e = (e). It was 
claimed experimentally that K ee (R) decays as a 

power law, 

K te {R)~R~», rj^R^L (1.5) 

with /r having a numerical value of 0.2 — 0.3. The an¬ 
alytic derivation of this law from the equations of huid 
mechanics, and the calculation of the numerical value of 


S n (R) ~ R Cn , 77 < R < L 


( 1 . 2 ) 
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the scaling exponent /.i have also been elusive goals of the 
theoretical research. 

In the present series of three papers we will present 
a systematic theory that culminates in the analytic esti¬ 
mate of the scaling exponents £„ and /x, in addition to 
the calculation of other quantities of experimental and 
theoretical interest. Since the road to these results is 
not straightforward, we present here a summary of the 
three papers in the series, hoping to motivate the reader 
to examine the necessary steps of analysis even though 
they are somewhat tedious. We believe that the proposed 
structure of the resulting scaling theory is rather beauti¬ 
ful and simple, justifying in part the lengthy path taken 
to reach it. 

The present paper which is numbered I deals with 
the perturbative theory of the correlation, response and 
structure functions of hydrodynamic turbulence. The 
main result of this paper (and see also Q) is that af¬ 
ter appropriate resummations and renormalizations the 
perturbative theory for these quantities is finite order by 
order in the sense that all the integrals appearing in the 
theory are convergent in the ultraviolet and the infrared 
limits . The meaning of this result is that there is no 
perturbative mechanism to introduce a length scale into 
the theory of the structure functions. When the per¬ 
turbative series diverges either in the infrared or in the 
ultraviolet the cutoff lengths that are used to tame the 
divergence serve as renormalization scales which can be 
used to change the scaling exponents from the prediction 
of dimensional analysis. This result indicates that as far 
as the perturbative theory is concerned there is no mech¬ 
anism to correct the exponents ( n from their K41 values. 
Of course, nonperturbative effects can very well furnish 
such a mechanism and therefore in paper II jf| (subtitled 
”A Ladder to Anomalous Scaling) we turn to the anal¬ 
ysis of nonperturbative effects. There is more than one 
source of nonperturbative effects that need to be stud¬ 
ied. One is obtained from summing up infinite perturba¬ 
tive series into analytic equations that display solutions 
that are not obtainable from the initial perturbation se¬ 
ries. Another one stems from sum rules or constraints 
on the perturbation theory as a whole. In paper II we 
study infinite resummations of ladder type diagrams. It 
is shown there that the renormalized perturbation theory 
for correlation functions that include velocity derivatives 
(to second or higher power) exhibit in their perturbation 
expansion a logarithmic dependence on the viscous scale 

| l],||. In this way the inner scale of turbulence appears 
explicitly in the analytic theory. The perturbative series 
can be resummed to obtain integrodifferential equations 
for some many-point objects of the theory. These equa¬ 
tions have also non-perturbative scale-invariant solutions 
that can be represented as power laws of r/ to some ex¬ 
ponents A. For example it will be shown in paper III 
H subtitled (Scenarios for Anomalous Scaling and In- 
termittency) that the correlation function of the energy 
dissipation field has such dependence: 


K ee (R) ~ (1.6) 

where A c = 2 — It has been argued || that if A < A c 
(a situation referred to as the “subcritical scenario”), 
then K41 is asymptotically exact for infinite Re. Then 
2 C 2 = (4 = 4/3 and the outer scale L disappears from 
( |l. 6 |) . In that case the exponent /x is identified with 
2(4/3 — A), and the renormalization length is the inner 
length 77 . In fact, in paper II it will be shown that the 
exponent A can be computed exactly, and that it takes 
on exactly the critical value A = A c . As a result of this 
the correlation K ee (R) can be shown to depend on R like 

K ee (R) ~ e 2 (i?/L ) C4 “ 2C2 . (1.7) 

In other words, the critical situation A = A c results 
in the disappearance of the inner renormalization scale 
and the appearance of the outer renormalization scale 
in ( 0 - In addition one notes that K ee (R ) decays as 
a function of R (i.e. the correlation is mixing) only if 
£4 < 2/2 which implies deviations from K41. Thus we 
will argue in paper II that the critical scenario A = A c 
goes hand in hand with multiscaling if K tt {R) is mixing 
and then /i is identified with £4 — 2/2 • Paper III which 
is the last paper in this series focuses on the analysis of 
the nonperturbative constraint furnished by the ballance 
equation which is obtained from the equations of motion 
of the structure functions. This analysis exposes a mech¬ 
anism for multiscaling which leads to the evaluation of 
the scaling exponents </ n from first principles. 


2. ON THE BACKGROUND OF THE ANALYTIC 
DESCRIPTION OF TURBULENCE 


Of foremost importance in the analytical approach 
have been those methods that rely on some per¬ 
turbative treatment of the Navier-Stokes equation 
Such treatments are fraught with diffi¬ 
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culties, which stem basically from two sources. The main 
source is the fundamental fact that turbulence has no 
small parameter. Thus, any naive perturbation scheme 
is bound to fail, being effectively a non-convergent series 
in powers of Re. Any truncation of such a series leads to 
disastrous nonsense. Evidently one must cleverly resum 
such series to obtain a formal perturbation expansion in 
terms of an effective coupling constant which differs from 
Re. The second source of difficulty is that the natural 
statistical objects which appear in the perturbation ex¬ 


pansions [hO 11 13J14] of the Navier-Stokes equation are 


correlation functions F a p{r, r', t, if) of the velocity field 
u(r, t) itself: 


T a p(v = (u a (r,i)up(r' ,t')) 


( 2 . 1 ) 


The problem is that the correlator F a p(r, r', t, if) is 
not universal, since it is dominated by contributions to 
u a (r,t) which come from the largest scales in the fluid 
flow which are determined by the features of the energy 


2 









injection mechanisms. This physical fact is reflected in 
the theory as infra-red divergences that have plagued the 
discussion of the analytic approach for decades. Indeed, 
all the early attempts to develop a consistent analytic 
approach to turbulence, notably the well known theories 
of Wylcl @ and of Martin, Siggia and Rose jn|], share 
this problem. 

It is ironic that the appropriate universal objects were 
known all along due to the insight of Richardson | [To|]T7| 
and later of Kolmogorov Jh| , who realized that velocity 
differences across a lengtliscale R do possess good prop¬ 
erties, since they are dominated by motions of scale R. 
In Kolmogorov’s language such a velocity difference is ex¬ 
pected to be universal when R is in “inertial range”. It 
seems advisable to develop the analytic techniques also 
in terms of similar universal objects. To this aim one 
may transform the equations of motion to new variables. 
A potentially promising approach could be to recast the 
Navier-Stokes equation in terms of Lagrangian coordi¬ 
nates as has been attempted by Kraichnan 0- The 
equation for the Lagrangian velocity involves however 
the Eulerian velocity which relates to the Lagrangian ve¬ 
locity via a complicated time evolution operator. A lot 
has been accomplished in this approach; firstly, Kraich¬ 
nan acheived a qualitative undertsanding of the different 
type of power law spectra that appear in different con¬ 
texts like scalar pasive diffusion, two and three dimen¬ 
sional turbulence, and Burgers dynamics. Secnodly, a 
numerical prediction for the Kolmogorov constant was 
obtained. However, the resulting perturbation theory 
does not have a diagrammatic representation in the sense 
that one does not know how to achieve infinite resumma¬ 
tions. The progress in understanding non-perturbative 
aspects remained unsatisfactory. 

Another approach that claimed success was due to 
Belinicher and L’vov Q who introduced an exact 
transformation to variables which they called “quasi- 
Lagrangian”. This term was not particularly successful; 
on the one hand it could give the impression of an approx¬ 
imate scheme, and on the other hand it implied proximity 
to the Lagrangian transformation. Both impressions are 
incorrect. The Belinicher-L’vov transformation is exact, 
and it has very little to do with the cumbersome features 
of the Lagrangian transformation. In terms of the Eu¬ 
lerian velocity u(r,f) they defined the held v(r 0 ,to|r,f) 
as 


v(r 0 , t 0 |r, t) = u[r + p(r 0 , f), f] (2.2) 


where 


p(r 0 ,t)=[ u[r 0 + p(r 0 ,r),r] . (2.3) 

Jto 

Note that p(ro,t) is precisely the Lagrangian trajec¬ 
tory of a fluid particle that is positioned at ro at time 
t = to. Nevertheless the held v(ro,fo|r,f) becomes in 
time very different from the Lagrangian held in which the 
Lagrangian trajectories of all space points are involved. 


The observation of Belinicher and L’vov was that the 
variables w(ro,fo! r >i) which are defined as 

w(r 0 , f 0 |r, t) = v(r 0 , f 0 |r, t) - v(r 0 , t 0 \r 0 , t) (2.4) 

exactly satisfy the Navier-Stokes equation in the limit of 
incompressible fluid. Consequently they could develop 
a diagrammatic perturbation theory in terms of these 
variables. The resulting theory was free of the two re¬ 
lated problems that we discussed above: the (ww) cor¬ 
relators were universal for |r — ro | in the inertial range, 
and the theory was free of infra-red (and ultra-violet) di¬ 
vergences. They also showed that the Kolmogorov 1941 
scaling is an order by order solution of the resulting the¬ 
ory. Two important properties of this formulation are: 
(i) The simultaneous correlators of v[ro|r,f] are identi¬ 
cal to the simultaneous correlators of u(r, t). The reason 
is that for stationary statistics the simultaneous correla¬ 
tors of an arbitrary number of factors of v[ro|r, t] does 
not depend on t, and in particular one can ta ke t = to- 
The property then follows directly from Eq.([2.3|). (ii) 
The correlators of w are closely related to the structure 
functions of u ( |l . l| ) . Clearly 

‘SnQr — r 0 |) = (|w(r 0 , f 0 |r, f)|") . (2.5) 

As the structure functions of the Eulerian velocity differ¬ 
ences were for years at the focus of experimental research, 
the formulation in terms of the variables w gives a direct 
link between theory and experiments. We will refer to 
the variable w(r 0 ,to| r 7 ^) as the BL velocity differences. 

For some reason these results have not been widely ac¬ 
cepted. One reason must have been the fact that the 
theory involved the subtraction of a velocity field at one 
special point (denoted ro above), and this, in the words 
of Kraichnan [jl9j , “cannot cure the infra red divergences 
that stem from far away regions”. Possibly the main 
reason for disbelief is that the Kolmogorov solution was 
found to be an exact order by order solution of the theory. 
There are two reasons not to believe a theory that repro¬ 
duces K41. Experimentally, it was widely stated that 
observations of the structure functions of velocity differ¬ 
ences are in contradiction with the K41 scaling, see for 
example 00 and references therein. Theoretically, it 
was believed since the days of Landau that the strong dis¬ 
sipation fluctuations “must” renormalize the exponents 
that are found within “naive” dimensional analysis (as 
is commonly found in the context of critical phenomena 
etc.). Lastly, the issue of the vertex renormalization was 
left out in 0 for future work. This also gave reasons to 
doubt whether the demonstration of the K41 scaling so¬ 
lution is to be accepted. The wisdom of the community 
about this issue was summarized recently in one sentence 
p2| : “The quasi- Lagrangian formulation serves only to 
provide a sophisticated diagrammatic reformulation of 
the Kolmogorov argument”. 

The aim of this paper is to clarify and solidify the the¬ 
oretical points that are needed to construct a Galilean- 
invariant and divergence-free renormalized perturbation 
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scheme for turbulence, and to answer all the doubts de¬ 
scribed above. Our strategy is to develop the theory from 
the start entirely in r, t representation in terms of new 
propagators which are correlatio ns a nd Green’s functions 
of the BL velocity differences (2.4); these objects pos¬ 
sess the good properties of being universal and of hav¬ 
ing transparent physical significance. We find that this 
new approach enjoys the merit that each diagram and 
the elements of the diagrams have definite physical inter¬ 
pretation. The simultaneous correlation functions of our 
theory are the functions of the Eulerian velocity differ¬ 
ences. We shall see that the physical reason for the lack of 
infra-red divergences is the local character of the hydro- 
dynamic interactions. This character is reflected in the 
fast decay of the Green’s functions when one of its spa¬ 
tial arguments goes away from a “ball of locality” whose 
radius is computed in the theory. The intuitive idea of 
the locality of interaction in fluids in r space finds an 
immediate reflection in the structure of the theory. Only 
interactions inside the ball of locality turn out to be im¬ 
portant. 

As a result of this theory we will explain why the K41 
scaling is indeed a consistent order by order solution of 
the theory. Notwithstanding, we will argue in papers II 
and III of this series that nonperturbative effects may 
furnish a multiscaling solution which differs from K41. 
The essence of the theory is that such a solution cannot 
be seen in order by order expansions. 

The structure of this paper is as follows: In Section 2 
we present the equations of motion, introduce the appro¬ 
priate variables, and derive, via the path integral formu¬ 
lation, the Dyson-Wyld equations for the second order 
correlators and Green’s functions (i.e the propagators). 
In Section 3 we discuss the physical meaning of the prop¬ 
agators and of the mass operators in our formulation, 
study the properties of the propagators, and show how 
they depend on time and on the coordinates. The cru¬ 
cial property of the Green’s function in BL variables is 
that its bare value is given as the difference between two 
projection operators. This serves as an initial value for 
the equation of motion of our Green’s function, leading 
to solutions that differ fundamentally from the Green’s 
function in the Eulerian representation. The properties 
discovered in Section 3 allow us to present in Section 4 
the proof of locality. This is the central section of this pa¬ 
per; it is shown that all the diagrams that appear in the 
infinite series for the mass operators in the Dyson-Wyld 
equations are local. This means that they have no di¬ 
vergences either in the large space- time separation limit 
or in the small separation limit, and they all contribute 
mostly in the “ball of locality” that is defined in section 
3. As a result we will be able to show that the only 
consistent scaling solution (order by order) has the K41 
scaling exponents. Of course, this approach suffers from 
the usual problem that although the series converges or¬ 
der by order, it may diverge as a whole. This is the 
possible source of nonperturbative effects that are being 
discussed in papers II and III of this series. Section 5 is 


reserved for a summary and a discussion. 


3. EQUATIONS OF MOTION, VARIABLES AND 
RENORMALIZATIONS 


A. Equations of motion 

The analytic theory of turbulence is based on the 
Navier-Stokes equation for the Eulerian velocity field 
u(r,t). In the case of an incompressible fluid they read 

du/dt + (u • V)u — i/V 2 u + Vp = f , V • u = 0, 

(3.1) 

where v is the kinematic viscosity, p is the pressure, and 
f is some forcing which maintains the flow. Since we are 
interested in incompressible flows, we project the longi¬ 
tudinal components out of the equations of motion. This 

is done with the help of the projection operator P which 
is formally written as P = — V _2 V x Vx. The applica¬ 
tion of P to any given vector field a(r) is non local, and 
it has the form: 


[Pa(r)] a = / dr'P a p(r — v')ap(r'), 


(3.2) 


where P a p{^ — r') is the inverse Fourier transform in k 
of the tensor P a p(k): 

/ rl 3 k 

^3 exp[— z(r-r') -k]P a/3 (k), (3.3) 

P a p{ k) = 6 a p - -j^kakp . (3.4) 


The calculation of the Fourier integral (3.3) gives 

P a p{r - r') = 8 a p5(v - r') (3.5) 
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r f 13 
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(r* - r' a )(rp - r' 0 ) 


r'|5 


Applying P to Eq. (0 we find 

(d/dt — vV 2 )u + P(u ■ V)u = P f . (3.6) 

We now perform the Belinicher-L’vov change of variables 


2.3) together with f[r + p(r 0 , t 0 , t), t] = 0[r o , f 0 |r, t]. 


Using the Navier Stokes equation and the chain rule 
of differentiation we find the equation of motion for 
v[r 0 , t 0 |r, t]: 


d 

dt 


r [r 0 ,-fo|r,f] + P[w[r 0 ,f 0 |r, t] ■ V]w[r 0 ,t 0 |r, t] 

= i/[V 2 + V 2 ]w[r 0 ,t 0 k,t]+P 0 (3.7) 


where w(r 0 ,to|ut) is the BL velocity difference defined 


in (2.4) and Vq acts on the ro variable. We note that 
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the equation of motion is independent of to, and from 
now on we drop the argument to in v(ro|r, t). Our aim 
is to develop a t heor y in terms of the BL field w. It 
can be seen from (3.7) that w satisfies the Navier-Stokes 
equation (with a diffe ren t forcing). However we will use 
the mixed equation (3/7) since the change of variables 
from u to v is more direct. 


B. Generating functional and perturbation theory 

The development of perturbation theory can be done 
using a number of e quiv alent formalisms, such as the 
Wyld diagrammatics |{n| or the path integral formula¬ 
tion Both approaches have been widely 

reviewed, and can be shown to give identically the same 
diagrams Jl5[ . We follow here the approach of Ref. 
||13||14|JT5| which is based on the generating functional 
Z{ A,m) 

Z( A,m) = -^ J Du(i)Dp(i) exp {il 

+ J dx[u{x) ■ A(x) + p(a;) • m(a:)]} , (3.8) 

where x = (r ,t). The quantity I is referred to as the 
effective action. It is customary to divide I into two 
parts, one quadratic and the other triadic in the field p 
and u: 


I = h 


Milt : 


(3.9) 

(3.10) 


Io = J dx\p a - vp a S7 2 u a ] 

J p a (x)D a p(x - y)p/ 3 (y)dxdy , 

lint = J dxp a (x)up(x)S7yu a (x) . (3.11) 

Here D a p(x — y) is the correlation function of a Gaussian 
random force f. We will not work in th ese varia bles, but 
rather make a chang e of variables in (3.8-3.11) accord¬ 
ing to ( 2.2-2~H) and (3.7). We will denote the resulting 
function as Z{ 1, m): 

Z( 1, m) = / DvTp exp {il 

+ J dx v(r 0 |a:) • l(:r) + p(a;) ■ m(a;)]} . (3.12) 

In the new variables the effective action can be again 
divided as in (13.91), but with 


I 0 = J dx\p a ^- - VPa^Wo] 

/ Pa ^' Da ^ x ~ y)pp(y) dxd V’ ( 3 - 13 ) 

dint ~ J dxp a (x)Wf3(x)\7pW a (x) . (3.14) 


Here ' D n n(x — y ) is the correlation function of the random 
force <fi in Eq.( |3.7| ) which is again assumed to be Gaus¬ 
sian (note that the initial force f is not Gaussian now, but 
this is of no consequence). In this formalism f DvDp is a 
functional integral over all possible realizations of the in¬ 
compressible velocity field v and the incompressible aux¬ 
iliary vector field p at all points of space and time. Note 
that in this formalism the auxiliary field p multiplies the 
equation of motion for v and that because of the restric¬ 
tion of incompressibility we do not have to display the 
transvers e pro jector in the effective action. The quantity 
Z in Eq.(3.12) is a normalization constant, 


Z= DvDp exp{i/} 


(3.15) 


which is equal to 1. The proof of this (2C]] is based on the 
expansion of the exp{*/}. The first term of this expan¬ 
sion is 1. The diagrammatic representation of the next 
terms will have closed loops made of Green’s functions; 
all these diagrams are zero because of c ausali ty. Upon 
expanding in powers of I m t we get from ( 3.12 ) a theory 
that contains summations of infinite terms compared to 
the theory with u(r, t ) as th e va ri able in the path inte¬ 
gral. To see this note that (2.2),(2.3) can be rewritten 
as 

u ( r , t) = v(r 0 |r — / v(r 0 |r 0 ,T)dr,t) . (3.16) 

Jt 0 

The RHS of this equation can be expanded in a Taylor 
series. The inversion of this series for v in terms of u is 
again an infinite series. We thus see that given powers 
of v contain infinite summations in terms of powers of u, 
and vice versa. It is extremely difficult to find the explicit 
diagrammatic resummation of the series in u to get the 
series in v. It is a great advantage of the path integral 
formulation that such infinite resummations can be ob¬ 
tained by a simple change of variables. This resummation 
is referred here as the Belinicher-L’vov renormalization. 
We reiterate that this renormalization is an exact formal 
step, and its usefulness needs to be demonstrated. It is 
alread y ob vious however that although the original func¬ 
tional (3F) produces c orrela tors which are not Galilean 
invariant, the new one ( |3.12 ) is free of this problem. This 
guarantees that the theory that is derived on the basis of 
this functional is Galilean invariant. 

An average of any product of the fields w and p is com¬ 
puted as a functional derivative of Z( 1, m) with respect 
to 1 and m in the appropriate order of differentiation, 
finally taking the limit l, m —» 0. For example 

F a p(r 0 \r,r',t,t') = (io Q (r 0 |r, t)wp(r 0 \r',t')) (3.17) 

8 2 Z( 1, m) 




\l,m —>0 


The perturbation theory is obtained by expanding ( 3. 12| ) 
in /; nt which contains the fields w and p. Since Iq is a 
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quadratic form in the fields, the natural objects which 
appear in the diagrammatic expansions are the correla¬ 
tors (ww), (wp) and (pp). The correlator (pp) is zero 
identically in this approach , whereas (wp) is 

the Green’s function, or the response to an additional 
external noise h(r, t) 

G a p(r 0 \r,r' (3.18) 
= i {w a {r o \r,t)p 0 {r ,f )} = ^ \ h ^ 0 . 

The first step in setting up the theory is the calculation 
of the bare Green’s functio n. Th e bare Green’s function 
G% is the response in w (|3.18| ) without the nonlinear¬ 
ity in the equation of motion. We can also introduce a 
Green’s function (7°^ which is the response of the field v 
without nonlinearity. From the equation of motion © 
we see that (7°^ satisfies the equation 

- uV 2 )G%(ry,t,t') = P aj3 {r - r')S(t - t'). (3.19) 


where the RHS stems from (3.2) with a(r) = ei5(r — 
r'), e is some unit vector. It is clear that (7°^ can be 
understood as the difference 

G a/ 3 ( r o|r,r',i,t') = G° a p(r,r',t,t') - G° af} (r 0 , r', t, t') . 

(3.20) 


From Eqs.( [3.1f - 3.20 ) we conclude that for t > t’ the 
Green’s function satisfies the equation 


J^-KVo + V 2 ) G° a0 (r o \r,r',t,t') 
= G° a/3 {r o \r 1 r , 1 0 + )S{t - t') 


(3.21) 


where 


G °/3( r ol r >r',0+) = P a p{ r - r') - P a p{ r 0 - r') . (3.22) 
The notation t = 0 + means the limit t —> 0 from above. 


The functions P a p (r—r') are given explicitly in Eq. (3.5). 
We recall that all Green’s functions vanish for time t < t 1 
due to causality. Note that the Green’s function at zero 
time is not proportional to a delta function in r — r'. 
This is a consequence of the condition of incompressibil¬ 
ity, which means that the speed of sound is infinite, al¬ 
lowing information to propagate in zero time. To derive 
the equation for the dressed Green’s function we note 
that the nonlinear part of the effective action has exactly 
the form of the Navier-Stokes nonlinearity, but in terms 
of v rather than w. Therefore the perturbation the¬ 
ory that this formalism produces is identical, in terms 
of graphical notation and in integral representation, to 
the standard perturbation theory that was developed in 
Refs. | |Toj . |TT}[i~3l|Tl . The important difference is that the 
objects that appear in the theory are correlators of w 
and responses of w, and the initial condition G°^(ro|r, 


r',0 + ) displayed in Eq. (3.23) is totally different from 
the initial condition in the Eulerian theory in terms of u. 
In the Eulerian theory G%(0 + ) has only one projection 
operator P a p{r — r') in contrast to ( 3.22| ). This difference 
will have far reaching consequences as we see later. At 
this point however we can use the standard procedure 
that results in the Dyson and Wyld equations for the 
line renormalized correlator and Green’s functions. The 
Dyson equation for the Green’s function reads 

^ - i/(Vg + V 2 ) G a p(r 0 \r,r',t) = G^rolr, r', 0 + )6(t) 

+ J dr 2 G° 5 (ro|r,r 2 ,0 + ) J dr x J dp 


x E, 5 7 (r 0 |r 2 ,iq, t 1 )G 7 ^(r 0 |r 1 ,r', t-t x ), 


(3.23) 


where t' was taken to be zero. Note that an immedi¬ 
ate consequence of Eq. ( 3.23| ) is that at short times 
lim t _» 0 + [G Q| g(r 0 |r, r', t)] = G^(r 0 |r, r', 0+). 

The Wyld equation for the correlation function has the 
form 

F a/3 (r 0 |r, r',i) = J dr x dr 2 J dMt 2 G a7 (r 0 |r, n, t x ) 

x [D lS (v 1 — r 2 , t — G + t 2 ) + d > 7 ^(r 0 |r 1 , r 2 ,t — t x + t 2 )] 

x Gi ( g(r 0 |r , ,r 2 ,t 2 ) • (3.24) 


In equation ( 3.23| ) the mass ope rator E is related to the 
“eddy viscosity” whereas in Eq.( 3.24 ) the mass operator 
$ is the renormalized “nonlinear” noise which arises due 
to turbulent excitations. The diagrammatic representa¬ 
tion for E and $ is well known, and is reproduced in 
Fig-1. 


4. PROPERTIES OF THE DRESSED 
PROPAGATORS 

In this section we present those asymptotic properties 
of the propagators which are used in the proof of locality 
in Sect.4. The reader who is not interested in the de tails 
of th e derivation ca n jus t pause to examine Eqs. (4.6- 
4.11) and then Eq. (4.51) before reading Sect.4. 


A. Rules for reading diagrams, and their physical 
meaning 

The objects in the diag rams of Fig.^j are the following: 
(i) the double correlator (3.16|) of the BL-velocity differ¬ 


ences w(r 0 |r, t), which is denoted by a wavy line that 
carries the designation for the points r x ,t x and r 2 , t 2 , see 
Fig-0- (ii) The Green’s function ( |3.18| ) which is repre¬ 
sented half by a straight line that denotes the field p, 
and half by a wavy line that stands for w. It also carries 
two space-time points, see Fig.0. Note that in contrast 
to the standard diagrammatic approach in the terms of 
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the Eulerian velocity the objects also carry a designation 
for the special reference point ro. (iii) The vertex T a p 1 
is: 

Eck/ 3 7 — ^07^/3 T ^a/3^7 • ( 4 . 1 ) 

In diagrams this vertex is denoted by a bold dot that is 
the junction of one straight and two wavy lines. This 
vertex is identical to the one which appears in the stan¬ 
dard theory in Eulerian variables (TyJTlJ] , and it is local 
in r, t space. The three lines that join at the vertex carry 
at that point the same space-time coordinate designa¬ 
tion. The differential operators in T Q ^ 7 act only on the 
propagators that meet the vertex with wavy lines. Every 
diagram appearing in the expansion for E Q , ( g(ro|r,tr', t') 
has an entry with a straight tail associated with r, t and 
an exit with a wavy tail associated with r ' ,t'. Dia¬ 
grams for < f> 0!/ 3 (ro|r, f, r', t') have two entries associated 
with straight tails carrying the designation r, t and r' 
,0 respectively, see Figj^. Every other space-time point 
that appears in the inner vertices of the diagrams is in¬ 
tegrated over. For example the one-loop diagrams for E 
and $ read: 


by the correlators of velocity differences between these 
points; if these correlators are relatively small, the con¬ 
tribution of the Green’s function to the total response is 
also small. 

The intuitive understanding of the Wyld equation is 
also straightforward. From the equation of motion CD, 
written schematically as 

W= [3/3t-^] { ~ |wV|w + ^' (44) 

we see that the nonlinear term [w • V]w is added 
to 0 and can be understood as an “additional non¬ 
linear noise” in the equation of motion. It is natu¬ 
ral to expect that the double correlator of the forc¬ 
ing will have a nonlinear contribution of the type 

([w(r 0 |ri, ti)Vi]w(ro|ri, fi)[w(r 0 |r 2 , t 2 )V 2 ]w(ro|r2,t2))- 

Indeed, the mass operator <f> Q , ( g(ro|r, r', t — t') can be writ¬ 
ten exactly as 

< f ) Q/ 3 (ro|r, r', t — t') (4.5) 

d d / i 

= w-o—TW*5( r o|ri,ri,r 2 ,r 2 ,fi,ti,f 2 ,f 2 ), 

or i 7 or2S 


d d 

^a/3(ro r, r ,t t ) (S ae b0*7 — )G ep (ro|r, r ,f f ) 


d d 

X (tpp-rp- + Sp6-pr) F 7 <5( r o|r, r', t-t'), (4.2) 

$ a/3 (r 0 1 r, r' ,t - 1') = ^ (S ae + 8 ai ) (6p p ^ 

d 


where F a ^p$ is the fourth order correlation of w. Thus 
Eq.( 3.24|) can be understood as the result of squaring 
Eq. ( I.4|) and averaging, up to the dressing of the Green’s 
functions. The content of the Wyld diagrammatics is 
the dr essin g of the bare Green’s function that appears 
in Eq.(T4), and the representation of the fourth order 
Fa-yps in terms of second order F a p. The dressing of the 
Green’s function results from the cross terms between 
+^p6-p-)Fjs{^o\ y, r ,t — t )F ep (ro\r,r ,t — t) . (4.3) [w • V}w and 0 in (4.4). Note that the analytic form of 


Note that the RHS of ([h^) is still an operator, since 
the operator d/dr'p acts on the Green’s function that ap¬ 
pears to its right in the Dyson equation ( |3.23 ). In con¬ 
trast, d ) a/3(ro| r , is a function, since the operators 

act only on the correlators that we see in Eq. (|4.3|) . At the 
one- loop o rder there are no inner vertices, and therefore 
(4.2) and (T3) do not have integrations. In general ev¬ 
ery diagram will have integrations over all intermediate 
vertices. 

The physical significance of the Dyson equation is that 
the dressed response is determined by hydrodynamic in¬ 
teractions involving intermediate points. For example, 
the response to forcing at r' has one direct, zero time 
contribution at r. However, for any finite time the re¬ 
sponse to forcing at r / is mediated by interactions at 
points ri which via E appear at point r. In the one- 
loop approximation E itself has a Green’s function that 
mediates directly between iq and r. In higher order con¬ 
tributions to E there are sequential contributions due to 
forcing at r' that are mediated by responses at various 
ri,r 2 ... until r is reached. E represents the dressed re¬ 
sponse which is the sum of all these sequential responses 
in the intermediate multiple sets of points. The Green’s 
functions that mediate intermediate points are weighted 


the $ a p in the one-loop order, Eq. (|4.3| ) follows directly 
from the Gaussian decomposition of F ai p$. 


B. Zero time properties of the renormalized 
propagators 

In the analysis of the asymptotic properties of the 
propagators it is convenient to translate the coordinate 
system such that ro = 0. We will denote the various 
quantities as G a( a(0|r, r', t), 3> a/ g(0|r, r', t — t') etc. The 
asymptotic properties of the bare Green’s function fol¬ 
low from Eq.( 3.2^ ), with the transverse projector given 
by (T5). Note that Gj^Olr,r',0 + ) is not a symmet¬ 
ric function of r and r' (also the finite time Green’s 
function is not symmetric!). The physical meaning of 
Ga/ 3 ( 0 |r, r', 0 + ) is the zero time difference between two 
responses of v at the points r and 0 to forcing at r'. For 
r' much smaller that r the main contribution has to come 
from the response at the point ro = 0: 


1 


0|r, r', 0+) « -P a/3 ( r') ^ — r' < r . (4.6) 

The next order contribution is about 1/r 3 which is much 
smaller. 
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In studying the other limit r <C H we note first 
that for r zero (i.e r = ro in arbitrary coordinates) 
^a/ 3 ( 0 | 0 , r', 0 + ) = 0 because of the exact cancellation 
of the two responses. For r small compared to r' we can 
expand the difference in responses in a Taylor series to 
find 


Gq^Olr, r ; , 0+) _ d 2 G a/3 (0\r,r\t) 


G° /3 (0|r,r',0 + ) 

^W Pap{jC,) 


r /5 


r«r 


(4.7) 


The last form demonstrates that in this limit 
r', 0 + ) is an odd function of r. 

The most important asymptotic properties of the cor¬ 
relator F a/ 3 ( 0 |r, r', t) for our discussion are those that ap¬ 
pear in its t=0 value. We can find these properties for 
time t = 0 by writing F a p(/d\r, r', t) identically as 

F a p( 0, r, r', 0) = S^r) + S a p( r') - S a p (r - r') (4.8) 
where 

S a p(r) = ([« a (r) - u a (0)][up(r) - 1 ^( 0 )]) . (4.9) 


The reason for ( |4.9| ) is that for t = 0 the correlators of v 
coincide with the correlators of u as has been explained 
in Sec. 2. Denoti ng S 2 (r) ~ S aa ( r) we find that in the 
regime rCr' Eq.(T8) reads 


F a /3(0\r,r', 0) ~ S 2 {r) +r 


dS 2 (r') 

dr 1 


(4.10) 


and 


Note that the function is symmetric in r, 
the opposit e lim it r /g> r’ we simply have to exchange 
r and r 1 in ( |f.lC| ). For t different from zero the correlator 
F a fj{0\r,r' ,t) is always smaller (by the Cauchy-Schwarz 
inequality) than its value at t = 0. This turns out to be 
all that we need for our discussion. The classical K41 


estimate of (4.10) is 

F a p(0\r, r',0) ~ e 


2/3 f 2/3 . r 


r n/3 


(4.11) 


where e is the mean energy dissipation per unit mass per 
unit time. 


C. Decay of the Green’s function in the short-time 
regime 

For the purpose of the proof of locality to be devel¬ 
oped in Section 4 we need to control the time integrals 
of the Green’s function. To this end we will need to 
estimate the characteristic decay time of G a /3(0|r, r', f) 
at long times. This is done in Sec. 3.D. Here we be¬ 
gin by estimating the characteristic decay time To ( 7 ', r’) 
of G Q ^(0|r,r',t) at short times t <C r 0 (r, r'). We will 
see that the first derivative of G a/ 3 ( 0 |r, r', t) vanishes at 
t = 0. Therefore we define ro(r,r') in the short time 
regime with the help of the second derivative 


f(r.r') 


dt 2 


|t=o+ 


(4.12) 


In principle ro(r, r') depends on the tensor indices a,/3. 
We will not be interested in numerical factors and we 
therefore disregard this dependence. From the point of 
view of the scaling theory of turbulence, typical time 
scales depend on the involved length scales. In K41 scal¬ 
ing the time scale r(r) is proportional to r 2 / 3 : 


r(r) 


gl/3 r 2/3 


(4.13) 


We are going to show that 

r 0 (r,r') ~ min{r(r), T(r')} = r(min{r, r 1 }) . (4.14) 

To compute the rate of chang e of G Q ^(0|r, r', t) we em¬ 
ploy the Dyson equation (|3.23). For length scales r and 


r i n the inertial range, we can neglect the viscous term 
in ( [l.23| ). The reason is that the viscous contribution to 
the time derivative ( |f,12| ) is of the order of i^/r 2 , which is 
much smaller than l/r(r). Initially we perfo rm th e eval¬ 
uation of the mass operator E on the RHS of ( 3.23| ) in the 
one-loop approximation, and then show that this evalu¬ 
ation is exact in the limit t —> 0. Substituting Eq.([h^) 
into ( |3.23 ) and writing the one-loop result schematically, 
(disregarding tensor indices) we find for t > 0 


d_ 

dt 


G(0|r, r',t) = J dr 2 G°(0|r,r 2 ,0 + ) J dri J dt\ (4.15) 

x— G(0|r 2 ,ri,ti) —[F(0|r 2 ,ri,ti)G(0|ri,r',t - G)]). 

Note that this schematic writing may be dangerous in 
that it obscures the possibility that some contributions 
might be odd under the operations rj —> —r! and/or 
r 2 —> — r 2 . To see this recall that the d/d r operations are 
odd whereas the rest of the integrand has even and odd 
contributions. Accordingly we will take the symmetry of 
the integrand into account when the result depends on 
it. 

The time arguments in F and G are bounded by t. 
For t <C r(r) we can neglect the time dependence of 
these propagators on the RHS of (4.15), and use their 
zero time values. Performing the (trivial) time integral 
and differentiating with respect to time we find 

J^G(0|r,r',t)«y dr 2 G°(0|r,r 2 ,0 + ) ( 4 - 16 ) 

xG°(0|r 2 , ri, 0 + )^-[E(0|r 2 , n, 0)G°(0|ri, r', 0 + )] . 

Next we note that the two-loop order has two more time 
integrations that for small t contr ibute to 9G(0|r, r', t)/dt 
starting at 0(t 3 ), and so to Eq.(lR) starting at 0(t 2 ). 
The n’th loop contributes starting from 0(t 2rl_1 ) to 
c>G(0|r, r', t)/dt. Thus we learn that for short times the 
expansion in the strength of interaction can be reformu¬ 
lated as an expansion in powers of the time. In the limit 























t —> 0 the one-loop approximation becomes exact. We 
can continue therefore to evaluate our quantities for short 
times in the one-loop order with impu nity. 

The analysis of the integral in ( |4.16| ) which we present 
next is quite cumbersome, due to the need to consider 
all the regimes of the relative sizes of r,r',ri and r 2 . In 
perform ing t his analysis we are going to use the K41 es¬ 
timate ( 4.11 ) of the correlator. We stress however that 
it will be shown that our conclusions are independent of 
the precise value of the exponents as long as they are not 
too far from the K41 exponents. This comment will be 
made quantitative later, when we understand the magni¬ 
tude of the window of locality that is def ined below. The 
result of our considerations is Eq. ( 4.36 ). 

Denote the minimum and the maximum of r, r' as 


G(0|r, r', t)\ t= o ~ r' dr 2 /rf /3 , wrong! (4.20) 

This integral indicates a divergence in the regime r\ ~ 
r 2 C A. This apparent divergence is not really there; it 
steins from disregarding the oddness of the integral with 
respect to a change in sign of iq and r 2 simultaneously. 
The odd sign ature that we missed is due to the r 2 • r'/r' 5 
[cf. Eq.(4/7)]. The rest of the integrand is even, and 
has been evaluated properly. Thus, looking at the next 
(even) order in the expansion of the last Green’s function 
we get instead of r 2 • r'/r 15 a term (r 2 • r') 2 /r . Finally 


instead of (4.20) we get the following convergent integral 


dt 2 


G(0|r,r',t)| t=o ~ r 


/—5 


dr 2 / r. 


1/3 
2 : 


right! (4.21) 


A = min(r, r'), A = max(r, r') . 


(4.17) 


I t wil l be shown that the main contribution to the integral 
( 4.16 ) comes from the regime r\ ~ r 2 ~ A. We begin by 
showing that the regimes in which n or r 2 or bo th ar e 
smaller than A contribute negligibly to the RHS of ( 4.16 ). 

In the regime iq < A < r 2 the relevant part of the 
integral over rq in (4.16) is 


The other cases, namely ?q <C r 2 < A and r 2 <C rq < 
A can be analyzed in the same manner with the same 
conclusion. In summary, we learn that the integrals over 
rq and r 2 in (4.16) may take A as their lower limit. 

Next we show that the upper limit of the integrals over 
iq and r 2 are bound ed by A. In the regi me ?q > A > r 2 
the relevant part of ( 4.16| ) is similar to (4.18), and may 
be evaluated as 


[ rfdriG°(0|r 2 , iq, 0 + )^-[F(0|r 2 , n, 0)G°(0|iq, r', 0 + )] 
Jo clr i 

(4.18) 


~*L drt 


1/3 ' 1/3 


L' 1 


r oo 

/ r?dnG°( 0|r 2 ,iq,0 + ) 

J A 

r) r 3 r°° 

x ^[ F (°l r 2i r ii°)G O (0|ri,r / ,0 + )] ~ 


(4.22) 

dr\ 


13/3 ' 


In the estimate of the LHS integral we used the evalu¬ 
ations G°(0|r 2 ,iq,0 + ) ~ 1/rf [cf. Eq. (|4.6|)1, <9/cbq 


To get the RHS we used the evaluations 
G°(0|r 2 ,iq,0 + ) - r 2 -iq/rf + (r 2 -ri) 2 /rf, F(0|r 2 , ri, 0) ~ 


1/n, F(0|r 2 ,iq,0) ~ [rf /3 + n/rq /3 1 [cf. Eq.(|4.11|)] and rf /3 + r 2 • iq/r 4/3 + (r 2 • n) 2 /rq 0/3 and G°(0|iq, F, 0+) 


G°(0|iq, i - ',) + ) ~ n/r' 4 [cf. Eq.(4/7)]. We took here the 


next order term in F because of the differentiation with 


respect to r 2 in ( 14.161) . W e see that in either case the 
main contribution to ( 4.18 ) comes from the upper limit 
rq ~ A. 

Second we co nside r the regime r 2 < A < r\. The 
relevant part of ( |4.16| ) now is 

[ r 2 dr 2 G°(0|r, r 2 ,0 + )^- [G°(0|r 2 , iq, 0 + )F(0|r 2 , n, 0)] 
Jo c,r ‘2 

(4.19) 


1/r' 3 . We need to consider higher order contributions 
since the leading order in F(0|r 2 , iq, 0) drops in the d/cbq 
differ entiation and the next order is odd. The integral in 
(4.22) contributes in the lower limit. 

In the regime r 2 > A > n the relevant integral looks 
like (4.19p, and is evaluated as 


dr 2 


r°° „ a 

/ r|dr 2 G°(0|r, r 2 , 0+) — [G°(0|r 2 , ri, 0 + )F(0|r 2 , ri, 0)] 
J A 

(4.23) 


dr 2 


11/3 


'A r-g 


4 and 


r 4 / 1/3 

'1 J 0 r 2 

where we used the evaluations G°(0|r, r 2 ,0 + ) 
1/rf, d/dr 2 ~ l/r 2 , G°(0|r 2 ,n,0+) ~ r 2 /r\ 
F(0|r 2 , ri, 0) ~ rf^ 3 . Again the contribution to the RHS 
integral is from the upper limit. 

The regime ri,r 2 < A is more delicate. In deed , 
the same kind of evaluation of the integral in ( 4.16| ) 
will go as follows: for rq ~ r 2 ), G°(0|r, r 2 , 0^ ~ 
1 /rf, G°(0|r 2 ,iq,0 + ) ~ 1/rf , F(0|r 2 ,ri,0) ~ r 2 /3 and 
the integral with respect to ri as rf. Using these evalu¬ 
ations we find 


In this case we evaluated G°(0|r,r 2 ,0 + ) ~ r • r 2 /r 2 + 
(r • r 2 ) 2 /rf, G°(0|r 2 , iq, 0+) ~ 1/rf, F(0|r 2 ,ri,0) ~ 
rf^ + ri •r 2 /rf /,3 + (ri •r 2 ) 2 /rf 0,/3 . We need the higher or¬ 
der terms for the same reasons as discussed above. This 
integral also contributes at A. 

Finally we need to consider the regime r\ ~ r 2 > A. 
The leading order evaluation is again zero by symmetry. 
Taking into account the next order contribution we eval¬ 
uate 


3 2 r 2 f c 

^G(0|r,r ,i)| t=0+ ~ ^ 


dr 2 

13/3 


(4.24) 
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In su mmar y, the contribution to the integrals over rq and 
r 2 in ( 4.1(j ) from regions much smaller than A and much 
larger than A are negligible. This is the first example of 
the property that we refer to in this pape r as “locality”. 
This property allows us to estimate ( 4.16 ) as the definite 
integral 


r(r, r') ~ (e) 1 / 3 r 2 / 3 ~ r(r), r<Cr'. (4.30) 


In the opposite limiting case, i.e. r r ', Eq.( 4.25| ) is 
employed with A = r' and A = r. After the evaluat ion o f 
the first and last Green’s functions the integral in ( 4.25| ) 
reads: 

d 2 


J^G(0|r,r',f)| t=0 + dr 2 r 2 G°(0\r, r 2 ,0 + ) j ^ 2 G(0|r, E, f)| t=0 +~ r' 3 ^ ^ (4.31) 

/ T O O 

dG?’?^G o ( 0 | r 2, ri ,0 + )—F(0| r 2, ri ,0) . 


x G>|r 2 , r 1; 0+) —[E(0|r 2 , ri , 0)G>|iq, r', 0+)]. (4.25) 


The integral (|4.25| ) offers us an easy evaluation of the 
characteristic time r(r, r') when r and r’ are of the same 
order, i.e. when A = 0(A). In this case we simply power 
count and find that 


d 2 1 

-^G(0|r,r ,f)| t=0 + ~ ^ 13 / 3 


Comparing this result with Eq.(4.12) shows that 


ro(r, r)~e 1 / 3 r 2 / 3 , 


(4.26) 


(4.27) 


which is identical to the Kolmogorov time r(r) given by 
( 4.13 ). We note that this finding is far from trivial; it 
requires the proof of locality. Once locality is shown, 
the analytical result must agree with dimensional anal¬ 
ysis as we indee d find . It is important to understand 
however that Eq.(4.25) offers more than can be achieved 
by dimensional analysis even with locality. For example 
we can discuss now tq(t, r') when r 3> r' and r <C r'. 
Since the ratio r/r' is dimensionless, dimensional anal¬ 
ysis cannot be employed to estimate the characteristic 
time r(r, r') in such cases. 

The f irst case to be discussed is r <C r'. We 
use Eq. (EH with A = r and A = r'. The first 
and last Green’s functions in ( |4.25| ) are evaluated as 
G°(0|r, r 2 , 0+) ~ r/rf and G°(0|iq, r', 0 + ) ~ ri/r' 4 . This 
leads to the integral 


dt 2 


G(0|r,r , ,t)| t=o + 


(4.28) 


r r rlrn F r) 

~-*Jr^lJr drir ^ G (0|r2 ’ ri ’0 +)F (0| r 2’ r 1 ’0)- 

We break now the integration over rq into the two do¬ 
mains (r, r 2 ) and (r 2 ,r / ). In the first domain the rq inte¬ 
grand is evaluated as 1 /r^ 3 , and as rq^/r 2 in the second 
domain. In both domains therefore the main contribu¬ 
tion comes for the regime rq ~ r 2 . Thus the rq integral 
is evaluated as 1 /ry . With this we can see that the r 2 
integral in ( |4.28| ) contributes mainly at r 2 ~ r, and can 
be evaluated as 1/r 4 / 3 . Finally 

J^G(0|r,r',f)[ t=0+ - ^ 73 , r < / . (4.29) 

Remembering that in t his regime G°(0|r, r', 0 + ) ~ r/r /4 
we conclude with (4.12) that 


We split again the ?q integral into two domains of integra¬ 
tion (r',r 2 ) and (r 2 ,r). In the first domain we evaluate 
G°(0|r 2 ,ri,0 + ) ~ 1/rf, and 9 2 [E( 0 |r 2 ,ri, 0)}/dridr 2 ~ 
l/r^ 3 . Substituting in the ?q integral we evaluate it as 
ln(r 2 /r , )/? , 2 ^ 3 . Evaluating then the r 2 integral we find 
that it contributes at r 2 ~ r'. We end up with 


J~2 G(CI|r, r ', f)| t=0 + 


r /13/3 


(4.32) 


In the second domain with r\ > r 2 we have to take into 
account two terms in the evaluation of G°(0|r 2 , iq, 0 + ), 
since the contribution of the first one will vanish by sym¬ 
metry G°(0|r 2 , iq, 0 + ) ~ r 2 • ri/rf + (r 2 ■ ri ) 2 /r\. Eval¬ 
uating 9E(0|r 2 , ri, 0)/9ri ~ r 2 /rj t;/3 we end up with the 
evaluation of the 7 q integrand as r 2 /rj 3 ^ 3 . The 7 q in- 
tegral contributes at ?q ~ r 2 and is evaluated as l/r 2 / . 
The r 2 integral contributes at r', and can be evaluated as 
1 /r' 4 / 3 . This means that the two domains of integration 
lead to the same evaluation ( 4.32| ) which is therefore the 
final evaluation in this regime of r > r'. Remem bering 
that here G°(0|r,r',0 + ) ~ 1/r' 3 we end up, using ( 4.12 ) 
with 


ro(r,r') ~ (e) 


^-i/ 3 r '2/3 


r(r'),r' < r . 


(4.33) 


The conclusion from Eqs.( 4.3C| ) and ( 4.33 ) is that in the 
short time regime the typical decay time of the Green’s 
function r(r, r') is always of the order of the smaller turn 
over times, as written in (4.14). This conclusion is intu¬ 
itively quite reasonable. 

It is very important to notice that we have prove n her e 
locality in a stronger sense than the one invoked in ( 4.25 ). 
It was shown that both the rq and the r 2 integrals take 
their main contribution from the regime rq ~ r 2 ~ A = 
min(r, r'). Although we used the Kolmogorov scaling for 
the evaluation of the integral, the result does not depend 
on that. We have a wide window of locality. We can 
change the scaling exponents of F( 0 |r 2 ,iq, 0 ) and/or of 
r(r) by a small number and still all our integrals will 
continue to be dominated by the same domains rq ~ 
r 2 ~ A. 


In fact, we can derive the result (4.14) for the typical 
time scale without invoking any numerical value for the 
scaling exponents. Assume that all the integrals in ( 4.25| ) 
are local in the sense that the largest contribution comes 
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from the region of integration n ~ r 2 ~ A = min(r,r'). 
Under this assumption we can evaluate G°(0|ri, r 2 ,0 + ) ~ 
1/A 3 , F(0 |r 2 , ri, 0) ~ S^/A) and d/dri ~ d/dr 2 ~ A -1 . 
Eq. (4.25) takes the form 


dt 2 


G(0|r, r', t) | t=0+ ~ AG°(0|r, A, 0+)S 2 (A)G°(0|A, r', 0+) 

(4.34) 


Using Eq. ( 1.1 2| ) we find 

G(0|r,r',0+) 


f(ur') 


AG°(0|r, A, 0+)S 2 (A)G°(0|A, r', 0+) 


. (4.35) 


Consider now the limiting case r <C r'. Now A = r 
and A = r\ and we can estimate G(0|r,r',0 + ) ~ 
A/A 4 , G°(0|r, A, 0+) ~ 1/A 3 and G°(0|A,r', 0+) ~ A/A 4 . 
Finally we get again (4.14) but with a more general ex¬ 
pression for t(A): 


t'o(u^) = t(A) , 
r(A) 




(4.36) 

(4.37) 


In the opposite case in which A = r r' = A all the 
Green ’s funct ions are evaluated as 1/A 3 and we regain 
( 4.36 ). Thus ( 4.36 ) is the final result for the characteristic 
time scale. This time scale is recognized as the typical 
turnover time of eddies of size A, having rms velocity of 
\J 52 (A). Clearly for K41 scaling we recapture the scaling 
law r(A) ~ A 2 / 3 , but in general our result is that the 
decay time of the Green’s function G(0|r, r'. t) is the turn 
over time of eddies of scale A = min(r, r'), independently 
of the exponent. Eq.(4.37) is one of the scaling relations 
between exponents of static and dynamic quantities that 
are used later to find the scaling solution of our theory. 

Note that the typical time scale r(r) becomes the only 
relevant time scale of G(0|r, r’,t) when r is about r'. In 
fact, when r ~ r' this is also the typical time scale of the 
correlator F(0|r,r',t), of the mass operator E(0|r, r',i), 
etc. The reason is that the problem contains no dimen¬ 
sional parameters, and therefore all the objects of the 
theory with the same spatial scale r must have the same 
characteristic time r(r). This is not an assumption, but a 
consequence of the structure of the theory and the prop¬ 
erty of locality that is proven in full in Sec.4. 


D. Characteristic decay time of the Green’s function 

In the previous section we evaluated the characteristic 
decay time To(r, r') of G(0|r, r', t) in the short time regime 
t <C To(r, r'). For the proof of locality we need to evaluate 
time integrals over G(0|r,r',f), and for this purpose we 
need to estimate the global decay time r(r, r') defined as 

/»00 

/ dt G(0|r, r', t) = r(r, r')G°(0|r, r', 0 + ) . (4.38) 

J 0 


Of course this definition makes sense only if the integral 
converges; we prove this fact next. For r and r' of dif¬ 
ferent orders we have no a priori reason to expect that 

r(r, r') is of the order of r 0 (r, r'). _ 

To prove the convergence of the integral in ( 4.38| ) we 
integrate Eq.(3.23) in time between t = — e and infinity. 
For r, r’ in the inertial interval we neglect the viscous 
term for the same reasons as discussed in Section 3C. We 
find 


r°° f) 

J gl G ^( °|r, r', t)dt = G^^ (0|r, r', 0+) 

+ J dt J dr 2 G° 5 (0|r, r 2 ,0 + ) J dr 1 J dt 1 

x E 57 (0|r 2 ,ri,fi)G 7/ 3(0|ri,r',f-fi) . 


(4.39) 


The Green’s function vanishes for t —> 00 , and it is iden¬ 
tically zero for all negative times. Accordingly, the LHS 
is zero, and after changing the order of integration over 
time on the RHS we arrive at 

G%{ 0|r, r', 0+) = - J dr 2 G° 5 (0|r, r 2 ,0+) Jdr 1 dt 1 

/»00 

X E 57 (0|r 2 ,r 1 ,t 1 ) dtG a/3 (0\r 1 , r', t) . (4.40) 

Jo 

Since the LHS is finite, both ti me in tegrals on the RHS 
converge. Thus the definition (4.38) is meaningful, and 
we proceed to estimate r(r, r’). 

To evaluate r(r, r') we forget the tensor indices, and 
then begin the analysis by assuming that the main con¬ 
tribution to ( |4.40p comes from the regi me r\ ~ r 2 . Eval¬ 
uating f dr 2 ~ r 3 , and making use of (4.38) we rewrite 


( pop 

G°(0|r,r , ,0+) = J rf dnG°(0|r, n, 0+) J 


dti 


xE(0|ri,ri,ti)r(ri,r , )G°(0|ri,r , ,0 + ) . (4.41) 


Since we proved that both time integrals in (4.40) are 
finite, the integral J 0 °° diiE(0|ri, ri, t\) can be evaluated 
as r(ri)E(0|ri,iq, t\ = 0). The time scale r(n) is the 
same as the one used in the context of the Green’s func¬ 
tion, because there is only one spatial scale available here. 
The mass operator E(0|ri,ri,ii = 0) can be evaluated 
in the one-loop app roxim ation ( |4.2| ) as S 2 {r)/r 5 . Using 
the scaling relation (4.37) this is equivalent to 


E(0|ri,ri,ti = 0) 


1 


rfr 2 (ri) 


Thus we can write 
G°(0|r,r , ,0+) 


(4.42) 


(4.43) 


r (n,r'), 


r 2 dr 1 G°(0|r,r 1 ,0+) V; G°(0|n, r', 0+) . 
nr 1 ) 


We are going to explore the implications of (4.43) on the 
allowed behavior of the ratio r(ri, r')/r(ri). As discussed 
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above, for r\ ~ r' this ratio should be of 0(1). For 
ri r' but both r and r' in the inertial interval we can 
expect a scaling behavior of the ratio 


- r(ri,r / ) 

r(n) 


r\ > r' 


(4.44) 


with some scaling exponent a. Similarly, for r\ <C r', we 
expect that 


T(ri,r') 

r(ri) 


P 


ri <C r' . 


(4.45) 


Consider now (4.42) in the domain (A,oo). The evalua¬ 
tion of the integrand is straightforward, and we write 


(4.46) 


[ x dr-[ 

V 

L r? 

n 


The meaning of this is deep: this equation says that 
the characteristic time of decay of the Green’s function 
G 7/ 3(0|r, r', t) is always bounded by the characteristic 
time of the first coordinate (r) which is the coordinate 
designating the difference between the two responses at 
r and 0. In the coordinate system with r 0 ^ 0 the decay 
of the Green’s function is bounded by the time scale of 
l r - r o|- 

We remember that up to now we considered the local 
contribution r\ ~ r^- There are good reasons to expect 
that this is the largest contribution. Howev er, if it turns 
out that this is not the case, then Eq.(4.40) implies that 
the estimate of the charact eristic time scale r(r, r') is 
even s maller than the bound ( 4.50| ) . Therefore the bound 
( 4.51 ) can be used safely. 


5. PROOF OF LOCALITY 


Convergence in the upper limit, which is guaranteed by 
the fact that the LHS is finite, requires a > — 1. In the 
domain (0, A) we evaluate 


G°(0|r,r',0 + ) 




(4.47) 


Convergence in the lower limit requires (3 < 1. With 
these restrictions (a > — 1 an d f3 < 1) the main con¬ 
tribution to the integral ( 4.43 ) comes from the domain 
( A, A) . For r » r' we evaluate the LHS and RHS of 
(4.43) as 


1 

r /3 


1 



(4.48) 


In order to satisfy this equation the integral must con¬ 
tribute in the lower limit. Thus a must satisfy a > 0. 
This is a stronger constraint than the one required by 
convergence. 

For r <?; r' we evaluate the LHS and RHS of (4.42) as 


r 



r dn 

V 

lr ri 

n 


(4.49) 


In this section we prove that all the diagrams for the 
Green’s function and correlators are dominated by inter¬ 
actions within a shell of a ball in physical space whose 
radius r will be determined later. The logic of the argu¬ 
ments is as follows. It is easy to evaluate an arbitrary 
diagram given that the coordinates of all vertices lie 
within a ball of radius r. The aim is to show that any 
contributions in which any of these vertices lies outside 
the ball of radius r is smaller that this “local” contri¬ 
butions. We will begin with the case in which only one 
vertex lies outside, and analyze the different geometries 
in which this may occur. We than consider the case of a 
group of vertices which lie close to each other but outside 
the ball, and finally the case where there are two sepa¬ 
rated groups of vertices lying outside the ball. This can 
be than generalized to any number of groups. The fact 
that all our integrals converge in both the upper and the 
lower limits serves as a proof that there are no perturba¬ 
tive mechanisms to renormalize the K41 exponents. This 
is the main point of this paper. 


A. The analysis of outer fragments of diagrams 


This integral must contribute in the upper limit in order 
to satisfy the equat ion. This re quire s f3 < 0. Combining 
these results with ( 4.44| ) and ( 4.45 ) we reach the final 
conclusion that 


r(r, r , )/r(7') < const (4.50) 

with const = 0(1) independently of the relation between 
r and r'. In terms of our goal which was the estimate of 
the integral over the Green’s function, we can write our 
final conclusion in the form 

7 /3(0|r, r', t) < r(r)G° /3 (0|r, r', 0 + ); always! 

(4.51) 



We begin by proving that if the coordinate of one ver¬ 
tex, say ri, in an arbitrary diagram is removed far away 
from all the rest of the vertices, the effect is to reduce 
the evaluation of the diagram as an inverse power of r\. 
Later we will use this result in our analysis of arbitrarily 
complicated diagrams. To facilitate the examination of 
the diagrams we will draw them in r, t space, such that 
vertices at a larger value of r will be positioned above 
vertices at smaller values of r, and larger times will be to 
the left of smaller times. Consider now a fragment of an 
arbitrary diagram in which one vertex whose coordinate 
ri is lying above all the other vertices in the diagram. 
In Fig.|] we show the three different possibilities of such 
a vertex at ri which is connected to another vertex at 
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smaller r via a Green’s function. The relative positions 
in space and time of the four adjacent vertices are clear 
from the pictures. We always take here r\ > r and our 
aim is to show that the major contribution to the inte¬ 
gral over ri comes from the region n ~ r. Note that 
due to causality the vertex at ri must be to the right 
of the vertex at r, and any other Green’s function that 
starts at ri has to point to the right. The wavy lines of 
the correlators can go either right or left, and we show 
only one representative example for each diagram that 
includes correlators. Since we are going to estimate the 
correlators by their t = 0 value, the time coordinates 
of the correlators do not matter. The three examples 
of Fig.|j are denoted as I a ,h and I c . The analytic ex¬ 
pressions for these three examples are shown below. For 
simplicity we do not display tensor indices: 




dti— G(0|r,ri,£i - 1) 
or 


(5.1) 


X 


— [F(0|ri,r 2 ,ti 
or i 


t 2 )F( 0|ri,r 3 ,fi -t 3 )] . 


We are going to bound I a from above. Accordingly, we 
take the largest possible value of the correlators, which 
is the t = 0 value . Aft er that we can bound the integral 
over t\ by using ( 1.51 ). We use the symbols ~ and oc in 
our evaluations, but we remember that we are actually 
estimating upper bounds; the quantities are smaller than 
what is displayed: 


r°° dri d 

I a (XT{r)J —[F(0|n,r 2 ,0)F(0|ri,r 3 ,0)] . (5.2) 


We now evaluate 

^-[F’(0|ri,r 2 ,0)F(0|ri,r 3 ,0)] « ^r^Tj/r^ 3 (5.3) 


where r :j is max{ r 2 , r 3 }, and = min{r 2 , r 3 }. Obviously, 
the integrals in (5.2) contributes at r 3 ~ r. 

The integral h, is 

h ~ J dri J diiJ^G(0|r,ri,fi - t) (5.4) 

d 

X —[G(0|ri,r 2 ,ti - t 2 )-F(0|ri,r 3 ,U -t 3 )] . 
or i 


In evaluating this integral we need to remember that in 
the generic diagrams the coordinate t 2 is integrated upon. 
The characteristic decay time of the Green’s function 
G(0|rir 2 ,fi — f 2 ) is bounded by r(ri). However because 
this Green’s function is connected to other propagators 
via vertex 2, the time restriction is at most r(r 2 ) and not 
r(r 3 ). Thus we write 


r°° dr^ d 

h oc t(t) / ——[G(0|r 1 ,r 2 ,0)F(0|r 1 ,r 3 ,0)] . (5.5) 

Jr r i or i 


We can evaluate now 


d 


-^[G(0|ri,r 2 ,0)F(0|ri,r 3 ,0)] « e 2/3 ^ . (5.6) 

1 r 3 r l 


We see that the integral (5.5) is local, contributing mostly 
at ~ r. 

Next we consider the integral I c : 


Ic ~ J dr! J dt 1 -^-G(0\r,r 1 ,t 1 -t) 
d 


(5.7) 


x — [G(0|ri,r 2 ,ti - t 2 )G(0|ri,r 3 , t x - t 3 )] 


where f; = min[ti,t 2 ]. In this case the integral 
/ G[0|ri,r 3 , (ti— t 3 )]dt\ does contribute a factor bounded 
by r(ri)G(0|ri, r 3 , 0 + ). We thus need to take into ac¬ 
count the factor r(r 3 ) in the integrand and write 


r (r) 


dn d 
r\ dri 


[r(n)G(0|ri, r 2 ,0)G(0|ri, r 3 ,0)] . 

(5.8) 


In the present case we use (4.6) and ij to write 


n • r. 


— [G(0|ri,r 2 ,0)G(0|ri,r 3 ,0)] ~ (5.9) 


i 1 1 


where r j is a unit vector in the direction of r j. Obviously 
also I c is local. 


B. Considerations involving arbitrary diagrams 


We can consider now arbitrary diagrams appearing in 
the Dyson equation. Every such diagram has a sequence 
of Green’s functions connecting an entry at r to an exit at 
r' via intermediate points r 3 , r 2 , ...r n , see Fig.|], panel (a). 
We call this sequence of Green’s functions the “principal 
path” of the diagram. From vertices along this principal 
path there emanate either a double correlator or a sub¬ 
tree that begins with a Green’s function that serves as 
its root. 

Consider first the situation in which r' < r. Our aim is 
to show that the main contribution to the diagram arises 
from the region of integration in which all the vertices are 
within the ball of locality of radius 0(r). Accordingly, 
we estimate first the mass operator E in the case that all 
the vertices are within the ball of locality. Using Eq. ( pd| ) 
the evaluation of E within the one-loop approximation is 
st r aight forward: 

E(r, r) ~ S 2 {r)/r 5 . (5.10) 


The two-loop diagrams contain two additional Green’s 
functions which contribute a factor of 1 /r 6 , one correlator 
that contributes a factor of S 2 (r) [cf.(2.5)], two vertices 
that contribute a factor of 1/r 2 , and two additional space- 
time integrations which give a factor of r 6 r(r) 2 . Together 
all these additional factors give S 2 (r)r(r) 2 /r 2 . Using the 
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scaling relation (4.23) we conclude that this is r°, and 
the two-loop diagrams are of the same order as the one- 
loop diagram. In fact, every order diagram in E gives 
the same evaluation as ( 0 - We remark that this is 
a common property of scale-invariant solutions in strong 
coupling systems. 

Next consider a contribution in which the coordinates 
of all the vertices in the diagram except one, say at a 
point r„, lie within a ball of radius r. This case can 
be dealt with using the analysis of sec.4.1: the situation 
locally must be one of the geometries (a), (b) or (c) in 
Fig.|. But we have shown that in all these cases the 
main contribution comes from the region r n ~ r. In 
fact the evaluation of the contribution will be diminished 
by a factor of (r/r„) a where the exponent a takes the 
following values according to the local geometry: 


f 7/3 (a) 
a = { 7/3 (b) 
{ 13/3 (c) 


(5.11) 


A situation in which a number of vertices go out of the 
ball of radius r but do not have direct in terac tions be¬ 
tween them is just a simple multiple of Eq.( |5.li ). There¬ 
fore the next situation to consider is that of a contribu¬ 
tion with an arbitrary number of vertices outside the ball 
of radius r with interactions. We take all of them to have 
r n of the same order, say of order r r. Such a blob 
of vertices can be connected to the ball of radius r via 
three legs again as in the situation (a), (b) and (c), or 
via a larger number of legs, see Fig.0. First consider the 
decoration of the uppermost vertex as in Fig.0 a,b and 
c. The decorated vertex has two additional Green’s func¬ 
tions, one additional correlator, and two space-time inte¬ 
grations. As we have learned such a combination scales 
like r°. Notwithstanding, the estimates of the reduc¬ 
tion in_the importance of th ese contributions compared 
to ( [5. 10 ) change from ( 5.11 ). The examples in Fig| a, 
b, c are evaluated as 


{S 2 (r)/r 5 }(r/r)‘ 2 


(5.12) 


The reasons for this evaluation are as follows: (i) The 
evaluation of the decorated vertex is 1 /r (recall that the 
bare vertex is d/df which has the same evaluation); (ii) 
the Green’s function connected via a straight line to the 
decorated vertex is evaluated as r/r 4 ; (iii) the space inte¬ 
gration gives f 3 . Lastly, the two additional l egs in Fig]|, 
a have no f dependence, and in toto we have ( [5.12|) . Note 
that in 5b and 5c we have one and two Green’s functions 
respectively. One could think that the time integrations 
over these Green’s functions would introduce factors r(f) 
and [i"(f )] 2 respectively. This is not so; the time integra¬ 
tion in the body of the diagram is restricted already by 
r(r) because of the propagators in the body. Therefore 
the appropriate factors to use are r(r) and [r(r )] 2 that 
do not change the dependence on r. 

Decorating the vertex further (which is equivalent to 
looking at higher order contributions to the renormalized 


vertex) d oes n ot change the estimate ( [5.12 ). The scaling 
relation ( 4.37 ) protects the evaluation of the contribu¬ 
tion. Therefore we consider now the situation in which 
the cluster of vertices at r contains an additional vertex 
that does not decorate the bare vertex, and which comes 
with an additional leg that connects the blob at r (Figjl| 
d and e). An additional vertex is worth 1 /r. The addi¬ 
tional Green’s function is worth 1/f 3 , and the additional 
space integration brings in f 3 . The additional time inte¬ 
gral is irrelevant since it is restricted by the shorter times 
in the body of the diagram. Adding more legs connected 
to the blob always reduces the contribution further. 

To complete the proof that the diagrams in the Dyson 
equation are dominated by their local contribution we 
now consider the situation that two blobs of vertices are 
located outside the ball of locality, say one blob at f\ 
and 7*2 7 see Fig!- There are two distinct possibilities: 
(i) the reduction factor depends on f 2 like some power 
of 7*2 and (ii) the reduction factor is independent of f 2 . 
In case (i) the power may be either positive or negative. 
If it is positive, then the main contribution comes from 
7*2 ~ fi, as shown in Fig. 0a. If the power is negative, 
the main contribution comes from f 2 ss r, as shown in 
Fig. 0c. Both situations have one blob outside the ball 
of locality, and are covered by the discussion above. In 
case (ii) we will gain a logarithm \og{ri/f 2 ) that will 
multiply the power (r/fi)“ that was derived above for 
the single blob situation. This of course does not change 
the conclusion that the main contribution comes from 
the ball of locality. 

The situation in which r and r' are of different orders 
of magnitude poses the question whether the ball of lo¬ 
cality is of radius r or r'. This question can be answered 
using the techniques developed above, but the natural 
framework to deal with this issue is to make use of the 
property of “rigidity” which will be defined in paper II of 
this series [0. The result || is that the ball of locality is 
always determined by the first coordinate in the case of 
the Green’s function G'(0|r, r', t) and by the smaller coor¬ 
dinate in the case of the correlation function F( 0 |r, r', t). 


C. Locality with respect to small intermediate 
coordinates and small separation distances 

So far we have only considered the possibility that in¬ 
termediate points Vj go out of the ball of locality, and we 
have shown that this diminishes drastically the contribu¬ 
tions to our diagrams. It is known that this is the danger¬ 
ous divergence that appears in the Eulerian perturbation 
methods. Nevertheless we discuss here for completeness 
the possibility that intermediate points r 7 approach ro, 
or the origin of coordinates in our system that is centered 
at r 0 . 

The discussion parallels to a large extent the steps that 
we took for the large distance case. Consider the situa¬ 
tions of Fig.0, but now since r\ is the smallest coordinate 
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in the group, one should invert the coordinate system. 
This is like looking at the figures upside down, and mir¬ 
rored across a line parallel to the r axis. (Causality al¬ 
ways forces G(0|r, n, ti — t ) to point from left to right in 
our coordinates where t increases to the left). The main 
difference with the previous calculation can be seen right 
away, say in the situation (a): 


6. DISCUSSION: THE IMPLICATIONS FOR 
SCALING SOLUTIONS 

We can examine now the implications of the fact that, 
order by order, our diagrammatic expansion is local. 
First we show that K41 scaling is indeed a solution. Tak¬ 
ing 


(a) = r(r) J r 2 ciri J^G(0|r, iq, 0+) (5.13) 

x A[ F ( 0 |r 1 ,r 2 ,0)F(0|r 1 ,r 3 ,0)] . 
or i 

Since the leading contribution to G(0|r, iq,0 + ), i.e. 1/rf, 
disappears with the d/dr di ffere ntiation, we take the next 
order which is 1/r 3 [cf. Eq. (|4.6| ) and the next line]. This 
factor does not influence the rq integration, and we gain 
a factor of rf compared with the integration in the do¬ 
main (r, oo). This difference now causes the integral to 
contribute at the upper limit, again at rq ss r. The same 
phenomenon recurs in all the relevant integrals, and we 
leave the completion of the proof, which is along the lines 
of the previous one, to the interested reader. 

Lastly, we need to consider also the situation in which 
two (or more) coordinates approach each other, e.g 
| r 7 ;,j = |iq — Tjl -C r. It is an immediate consequence 
of our choice of BL variables that nothing depends ex¬ 
plicitly on this separation. Since we integrate over this 
separation from zero to r, the phase volume rfjdrij en¬ 
sures that the main contribution comes from |iqj| ~ r. 

In summary, we have shown that the exact resumma¬ 
tion of the Eulerian perturbation theory which results 
from the Bclinicher-L’vov transformation leads to a new 
perturbation theory in which all the diagrams are local. 
We reiterate that “locality” here means that the major 
contribution comes from intermediate coordinates that 
are in a shell of radius r. The width of this shell is also 
of the order of r. This was proven in the last paragraph 
when we showed that the contributions from the center 
of the ball of locality with rq <C r are also negligible. If 
we measure separation distances on a logarithmic scale, 
as is appropriate in a scale invariant situation, we can 
state that the major contribution to the diagrams come 
from the narrow shell of the ball of locality. 

We can state these results in an intuitive fashion. With 
a perturbation at a point r', from observing the differ¬ 
ences of responses at points r and ro, we learn from our 
diagrammatic expansion that the coordinates of all the 
vertices [which designate the position of the hydrody¬ 
namic interactions (u • V)u] that appear in any order of 
the expansion must lie within a physical ball of fluid of 
radius |r — ro | in order to have a significant effect. Also 
the double correlator is affected by intermediate points 
in a finite ball of fluid of radius max(jr — Tq |, | r' — ro |) - 
This is an analytic demonstration of the locality of the 
interactions between turbulent fluctuations in physical 
space. 


S 2 {r) 


„(2 


r (r) 


( 6 . 1 ) 


we see that Eq.(4.37), which has been obtained solely 
from the property of locality, furnishes one scaling rela¬ 
tion between (/ 2 and z, i.e. 


C 2 + 2 z — 2 . 


( 6 . 2 ) 


We can easily find a second scaling relation by using our 
theory to calculate the 3-point correlation function 


F a p~/{ro\ri,r 2 ,r 3 ,ti,t 2 ,t 3 ) (6.3) 

= (w a (r 0 |ri,fi)wg(ro|r2,f 2 )«> 7 ( r o|r3,£3)) • 

The simultaneous correlation F , Q ,/3 7 (r 0 |r, r, r, t, t, t) is 
identical to the 3- rd o rder structure function of the Eule¬ 
rian field, cf. Eq.( [2.5|) , and it is well known that for small 
viscosity it scales (in the inertial range) like r. The 3- 
point correlation function has a totally analogous theory 
to the one described above for the 2 -point function, and 
a typical first order diagram for this quantity is shown 
in Fig.|[ By repeating the analysis that led to the proof 
of locality we see that the main contribution to the inte¬ 
gral f dridti Fig. 7a comes from the region rq ~ r and 
ti ~ r(r) . Consequently we can evaluate S 3 (r) in the 
first order as 


S 3 (r) ~ T{r)[S 2 (r)} 2 /r . 


(6.4) 


The next order contributions to F a p^ are nothing but 
the renormalization of the bare vertex, see for example 
Fig. @b. One can prove following the ideas of Sec. 4 that 
all these diagrams a re loc al. The property of locality and 
the scaling relation (4.37) or equivalently ( |6.2|) imply that 
all the n’th order diagra ms have the same evaluation as 
the first order diagram ( |6.4|) . Since S 3 (r) ~ r, we find a 
second scaling relation 


[S 2 (r)] 2 T{r) ~ r 2 . 

In terms of the scaling exponents this is 
2 C 2 + z = 2 . 


(6.5) 


( 6 . 6 ) 


The solution of (6.2) and (Eh5) is ( 2 = 2/3, z=2/3, which 
are the standard K41 exponents. 

The usual mechanisms to renormalize exponents are 
connected to the divergence of the integrals appearing 
in higher order diagrams. The cutoff lengths that are 
needed to control the divergence can be used to form 
dimensionless corrections to the scaling functions, and 
thus to renormalize the exponents. Since our diagrams 
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are local, there are no typical scales which appear in this 
way, and we do not have a simple way to renormalize 
the scaling exponents of the structure functions. Indeed, 
if there is a renormalization of the exponents, it must 
stem from nonperturbative effects. The study of these 
effects and the elucidation of nrultiscaling is the subject 
of papers II and III in this series. 
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APPENDIX A: RELATION BETWEEN THE 

CORRELATION FUNCTIONS BEFORE AND 
AFTER THE BL TRANSFORMATION 

In this Appendix (which follows ^8|) we examine the 
relation between the BL transformed correlation func¬ 
tions and the correlation functions of the Eulerian fields. 
The aim of this appendix is to clarify how the stationarity 
and homogeneity of the original fields are reflected in the 
properties of the transformed fields. We introduced the 
BL velocity v by the expression ( |2.2| ): v(f, r) = u(t, r+p) 
where u is the Eulerian velocity and p(f) introduced by 
( [2.3P is governed by the equation 

dp/dt = u(f,r 0 + p) = v(f,r 0 ), p(t = t 0 ) = 0. (Al) 

Here to is a marked time and r 0 is a mark ed point which 
can be chosen arbitrary. Equation © describes the 
Lagrangian trajectory ro + p of the liquid point starting 
from the point ro at f = to. 

Let us introduce a designation for the product of the 
BL velocities taken in different spatial-time points, 

^{p} = v(t 1 ,r 1 )v(t 2 ,r 2 ).. .v(t n ,r n ) (A2) 

= u [ti, ri + p(ti)] u [t 2 , ri + p(t 2 )] ■ • ■ u [t„, r n + p(t n )}. 

Next we transform the expression for the different-time 
correlation function of BL velocities, (T{p}), into the 
form which is convenient for utilizing the relation be¬ 
tween u and v. For this goal we will treat p as an i nde- 
pendent field taking into account the equation © via 


('F{p}) = J DpD/3 exp J dtfddp/dt 

x ^exp (-i J df/3u(f,r 0 + p)^j T{p}^> ■ (A4) 


In the above formal transformations we omitted the 
functional determinant arising at passing to integration 
over p, j3. To recognize the value of the determinant we 
shou ld first discretize the problem introducing instead of 
©) a discrete scheme, e.g. 


p{t n + 1 ) ^ p{t n ) - u (t, r 0 + p{t n )) At = 0 , 
t n =to+ nAt. (A5) 


Now the functional 5-function in (A3) is the product of 
usual 5-functions wit h arg uments defined by the left-hand 
sides of the equation (A5), the product being taken for all 
numbers n /0 (since in accordance with (Al) p(fo) = 0 ). 
The functional integration in (|A3|) is now the product of 
integrations over p(t n ), n / 0. We see that for n > 0 
the consequent integration over p{t n ) actually gives unity 
since 


j dp{t n+ x)8 (p(t„+i) - p(t n ) -...) = 1 • 

But when we go over n in the negative direction the in¬ 
tegrals over p{t n ) produce nontrivial factors associated 
with a dependence of u on p. The product of these fac¬ 
tors is just the functional determinant noted above. 

Usually this problem is avoided by putting to into — oo, 
after that the actual region of t produces the unitary de¬ 
terminant. It is possible to do only if we deal with a 
situation as a stationary one which is not sensitive to 
initial conditions. Just this reason explains that the cor¬ 
responding functional determinant can be regarded to be 
equal to unity in the path integral formulation of the 
Wyld-Martin-Rose-Siggia technique. But our problem 
is in proving time homogeneity of BL correlation func¬ 
tions, therefore to is a finite time and we should consider 
times t both before and after to- Fortunately in our case 
the incompressibility condition ensures the unitary of the 
determinant. The point is that for small At in the ap¬ 
proximation needed for us the integrals 

J dp(t n -i)S ( p(tn ) - p(tn-l) — •••) = 1 — AfVu 
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are equal to unity due to the incompressibility condition. 
Thus we see that the condition is of the crucial impor¬ 
tance for the Belinicher-L’vov scheme. 

Taking into account the definition (Al) we con¬ 
clude that the expression in the angle brackets in (A4) 
is an average of a combination of Eulerian velocities 
u [t, r + p(t)) with a function pit). This average is obvi¬ 
ously translationally invariant. It allows one to shift all 
space arguments at a constant vector R. L et us choose 
R = p(r) where r is some time. Then ( |A4| ) acquires the 
form 


(vE'Ip}) = j Vp'VPexp J dt/3dp'/dt 
x ^exp ^-i J dt(3u(t, r 0 + p')^ , (A6) 

where p'(t) = p(t ) — p(t). Note that p'{r ) = 0 . Integra¬ 
tion over p,(3 in (A6) gives 




(A7) 


where p'[t) is governed by the same equation & as p 
but with the ini tial condition p'(r) = 0. 

The relation (A7) means that BL correlation functions 
are stationary in time since th ey d o not actually depend 
on the marked time r. Using (A7) one can easily prove 
the coincidence of simultaneous correlation functions of 
BL and Eulerian velocities. Namely, taking i n (|A2|) all 
times t\, t 2 , • • ■ to be equal to t and choosing in (|A7p r = t 
we immediately find 


(v(t,ri)v(t,r 2 ) ■ ..v(t, r„)) 

= {u{t, v{)u{t, r 2 ) ... u(t, r„)), (A8) 

since p’(t) = 0 in this case. Although we examined corre¬ 
lation functions of BL velocities, all arguments presented 
above can be applied to arbitrary fields e.g. pumping 
forces. It implies e.g. that BL Green’s functions are 
stationary in time and that their simultaneous values co¬ 
incide with ones of Eulerian Green’s functions. 


* L’vov’s e-mail: fnlvov@wis.weizmann.ac.il. 

* Procaccia’s e-mail: cfprocac@weizmann.weizmann.ac.il. 

[1] A.S. Monin and A. M. Yaglom, Statistical Fluid Mechanics 
(MIT press, 1975) 

[2] K.R. Sreenivasan, Annu. Rev. Fluid Mech. 23,539 (1991). 

[3] M. Nelkin Advan. in Phys. 43 143-181 (1994). 

[4] V. I. Belinicher and V. S. L’vov. A scale-invariant 
theory of fully developed hydrodynamic turbulence. 
Sov. Phys. JETP, 66:303-313, 1987. 

[5] V. S. L’vov and I. Procaccia. Exact resummations in the 
theory of hydrodynamic turbulence. II. A ladder to anoma¬ 
lous scaling. Phys. Rev. E, 1995. Next paper in this issue. 


[6] V.S. L’vov and I. Procaccia, Exact resummations in 
the theory of hydrodynamic turbulence. III. Scenarios 
for Anomalous Scaling and Intermittency, Submitted to 
Phys.Rev.E 

[7] V. S. L’vov and V. V. Lebedev. Anomalous scaling in 
Kolmogorov 1941 turbulence. Europhys. Lett., 29 681-686 
(1995). 

[8] V. V. Lebedev and V. S. L’vov. Scaling of correlation func¬ 
tions of velocity gradients in hydrodynamic turbulence. 
JETP Letters, 59:577-583, 1994. 

[9] V. S. L’vov and I. Procacia. “Intermittency” in turbu¬ 
lence as intermediate asymptotics to Kolmogorov’41 scal¬ 
ing. Phys. Rev. Lett., 74, 2690 (1995). 

[10] H. W. Wyld. Formulation of the theory of turbulence in 
an incompressible fluid. Ann. Phys., 14:143-165, 1961. 

[11] P. C. Martin, E. D. Siggia, and H. A. Rose. Statistical 
dynamics of classical systems. Phys. Rev. A, 8(l):423-437, 
July 1973. 

[12] R. H. Kraichnan. Eulerian and Lagrangian renormalization 
in turbulence theory. J. Fluid Mech., 83(2):349-374, 1977. 

[13] C. DeDominics. J. Physique (Paris), 37:Cl-247, 1976. 

[14] H.K. Jansen. Z. Phys. B, 23:377, 1976. 

[15] V. S. L’vov and I. Procaccia. Exact resummations in 
the theory of hydrodynamic turbulence. 0 Line-resummed 
diagrammatic perturbation approach. In F. David and 
P. Ginsparg, editors, Lecture Notes of the Les Houches 
1994 Summer School, 1995. In press. 

[16] L. Richardson. Weather Prediction by Numeric Process. 
Cambridge LTniv. Press, Cambridge, 1922. 

[17] Lewis F. Richardson. Atmospheric diffusion shown on a 
distance neighbor graph. Proc. R. Soc. Lond. A, 110:709- 
737, 1926. 

[18] A. N. Kolmogorov. The local structure of turbulence in 
incompressible fluid for very large Reynolds numbers. 
C. R. (Doklady) Acad. Sc. URSS, 30:9-13, 1941. 

[19] R. H. Kraichnan, 1994. Private communication. 

[20] G. Stolovitzky and K.R. Sreenivasan. Scaling of structure 
functions. Phys. Rev. E, 48(1):R33-R36, 1993. 

[21] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, F. Mas- 
saioli, and S. Succi. Extended self-similarity in turbulent 
flows. Phys. Rev. E, 48(1):R29-R32, July 1993. 

[22] C-Y. Mou and P. B. Weichman. Multicomponent turbu¬ 
lence, the spherical limit, and non-Kolmogorov spectra. 
CalTex, , 1994. Preprint. 

[23] V. S. L’vov, I. Procaccia, and A. Fairhall. Anomalous 
scaling in fluid dynamics: The case of passive scalar. 
Phys. Rev. E, 50(6):4684-4704, December 1994. 

[24] V. S. L’vov and I. Procaccia. Correlator of velocity differ¬ 
ences and energy dissipation as element in the subcritical 
scenario for non-Kolmogorov scaling in turbulence. Euro¬ 
phys. Lett., 29 (4) 291-296, 1995. 

[25] V. S. L’vov and V. V. Lebedev. Exact relations in 
the theory of fully developed hydrodynamic turbulence. 
Phys. Rev. E, 47(4)4794-1802, 1993. 

[26] E. Kats and V. Lebedev. Fluctuation effect in the dynamic 
of liquid crystals. Springer Verlag, Heidelberg, 1993. 

[27] A. A. Praskovsky and S. Oncley. Correlators of velocity 
differences and energy dissipation at very high Reynolds 
numbers. Europhys. Lett., 29 (9) 635-640 1994. 

[28] V. V. Lebedev and V. S. L’vov, e-print in Los-Alamos e- 


17 









board No.9410003 (1994) chao-dyn@xyz.lanl.gov 


W — AAA 5 w 0 — aaa , 


FIG. 2. Diagrammatic representation of the renormalized 
series expansion for the mass operators, (a) The mass opera¬ 
tor S of the Dyson equation (2.23) and (b) the mass operator 
4? of the Wyld equation (2.24). 



F ~ AAAAAA, Fo ~ AAAAAA, 

G ~ Ann_, G 0 ~ nnn-. 

FIG. 1. Graphical notation for the perturbation expansion. 
The symbol used are the following: short wavy lines stand for 
the fluid velocity (think about waves). A thin wavy line stands 
for uo, whereas a bold wavy line represents the full solution 
u(r, t). A straight line stands for the field p(r, t ) that appear 
in the functional integral approach. The Green’s function, 
which is the response in the velocity to some force is made of 
a short wavy line and a short straight line representing the 
force. Again this symbol appears in thin and bold variants, 
the former stands for the bare Green’s function (2.20), and 
the bold for the dressed Green’s function (2.18) The vertex 
(3.1) is a fat dot with three tails. One straight tail belongs to 
the Green’s function, and two wavy tails stand for velocities. 
A long wavy line will represent correlation functions of veloc¬ 
ities. Again, the thin and bold variants are bare and dressed 
correlators respectively. 


(a) 




FIG. 3. Three different arrangements of four adjacent ver¬ 
tices in a fragment of a diagram, drawn in r, t space. 






(b) 





FIG. 4. Fragments of diagrams which are used in the proof 
of locality. 
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FIG. 5. Fragments of diagrams which are used in the global 
considerations in the proof of locality. 



FIG. 6. Diagrams in the expansion for the triple correlation 
function, (a) All the first order contributions, (b) some of the 
third order contributions. 
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